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Abstract 

Using first-principles calculations, we systematically study the atomic structures and electronic 
properties for two dimensional triangular ZnO nanoislands that are graphite-like with monolayer 
and bilayer thickness. We find that the monolayer ZnO nanoisland with O terminated zigzag edges 
is magnetic at its ground state, with the magnetism coming from the O edge states. The other 
monolayer and bilayer ZnO nanoislands with different edge structures are all nonmagnetic at their 
ground states. It is further revealed that for different ZnO nanoislands, their magnetic properties 
are quite dependent on their sizes, with larger nanoislands having larger magnetic moments. 
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As a semiconducting metal oxide with a direct band gap, zinc oxide (ZnO) has vast 
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tinctures of different morphologies have been fabricated and studied |5[. Especially, 
graphite-like hexagonal ZnO nanofilms have been successfully fabricated [6] and extensively 
studied for possible usages in electronic devices [3, S Is]. Moreover, similar to the situations 
in graphene nanostuctures, ZnO nanotubes and nanoribbons have then also been investi- 
gated fo.. possible potential appl,ca«o„.fl[l HQ. It is noted that the monolayer 



ZnO nanoribbons with zigzag edges exhibit magnetic behavior [13|. Besides carbon nan- 
otubes and graphene nanoribbons, two dimensional (2D) graphene nanoislands have also 
jeen paid vast attentions because of their unique flexibility in tuning magnetic properties 
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151]. Thus when turn to ZnO nanostructures, one may wonder whether 2D ZnO nanois- 



lands show tunable magnetic properties or not. Especially when ZnO nanoisland samples 
is experimentally attainable by depositing hexagonal ZnO nanofilms on the Ag(lll) surface 
61]. Motivated by this question, in this paper we present systematical investigations on the 
atomic structures and electronic properties of 2D triangular ZnO nanoislands. 

Moreover, the so-called (i°-magnetism have recently attracted extensive interests, which 
represents for magnetic behaviors in semiconducting materials (including ZnO) with the 
absence of transition metal doping. It was firstly observed in Hf02 thin films where defects 
or oxygen vacancies cause magnetic behaviors [l^. Later, it has also been theoretically 



predicted that carbon doping in ZnO [17|] and cation-vacancy in GaN and BN ISj both 
will lead to magnetic behaviors. In addition to these previous researches, here through 
first-principles calculations, we further reveal that the monolayer ZnO nanoisland with O 
terminated zigzag edge also shows magnetism in its ground state, even without any vacancies 
or defects. 

In the present work, the first-principles calculations are carried out using the DMOL^ 
package. Density functional semicore pseudopotentials (DSPP) [l9] are adopted to replace 
the effects of core electrons of Zn and O with a simple potential including some degree of rel- 



ativistic correction. The spin-polarized PW91 2(| functional based on generalized gradient 



approximation (GGA) is employed to take into account the exchange and correlation effects 
of electrons. For valence electrons, the double-numerical basis with polarized functions is 
adopted for expansion of the single-particle wave functions in Kohn-Sham equations. For 
better accuracy, the octupole expansion scheme is adopted for resolving the charge density 
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and Coulombic potential, and a fine grid is chosen for numerical integration. The charge 
density is converged to lxlO~^ a.u. in the selfconsistent calculation. In the optimization 
process, the energy, energy gradient, and atomic displacement are converged to 1x10^^, 
1x10^^ and 1x10^^ a.u., respectively. The atomic charge and magnetic moment are ob- 
tained by the Mulliken population analysis. 

In the experiment by Tusche et al, the observed triangular 2D ZnO nanoislands all have 
the size of about 20 A 6]. So we here start by investigating several kinds of ZnO nanoislands 
at this size, including the monolayer and bilayer ZnO nanoislands with both armchair and 
zigzag edges. It should be noted that there are two different kinds of zigzag edges (respec- 
tively the Zn- and 0-terminated) for monolayer ZnO nanoislands. In total, five different ZnO 
nanoislands are studied here: the monolayer armchair-edge ZnO nanoisland (AZnONI), and 
zigzag-edge ZnO nanoisland with Zn- (Zn-ZZnONI) and 0-terminated edges (0-ZZnONI), 
the bilayer armchair-edge (BAZnONI) and zigzag-edge ZnO nanoisland (BZZnONI). The 
optimized atomic structures for these ZnO nanoislands are all shown in Fig. 1. 

During geometry optimizations, different kinds of ZnO nanoislands show different struc- 
tural reconstructions. As shown in Figs. l(a)-(c), after geometry optimizations for the 
monolayer ZnO nanoislands, all the atoms are still inside a same plane. For the AZnONI, 
the atoms at the edges and corners are distorted from their equilibrium places, and nine new 
Zn-Zn bonds form at the three edges as shown in Fig. 1(a). For the monolayer zigzag-edge 
ZnO nanoislands, the structural reconstructions only occur at the three corner areas. One 
can see from Fig. 1(b) that in the Zn-ZZnONI structure, the two Zn atoms at each corner 
come close with each other and form a new Zn-Zn bond, with the three O atoms at the cor- 
ners pushed out from their original sites. In the 0-ZZnONI structure shown in Fig. 1(c), the 
two O atoms at each corner move away from each other a little making the three Zn atoms 
at the corners pulled in from their original sites. The reconstructions for the BAZnONI 
and BZZnONI structures are more complex. Figures. 1(e) and (g) are respectively the side 
view for the atomic structures in the blue and red squares in Figs. 1(d) and (f), from which 
one can see that the atoms at corners and edges are no longer in the same plane with the 
middle atoms. The bond angles at these distorted areas are decreased to be smaller than 
120°, indicating that these atoms incline to transit to the wurtzite configuration. 

The electronic properties for the 2D ZnO nanoislands are then studied, and found to 
depend critically on their edge structures. Figure 2 shows the calculated energy levels for 
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the five geometrically optimized ZnO nanoislands with different edge structures. From Figs. 
2(a), (d), and (e), we can clearly see that there exist large energy gaps for the AZnONI, 
BAZnONI and BZZnONI structures. This is in accordance with the saturation that the 
ratio between Zn and O atom numbers is 1:1 in them. The Zn-ZZnONI and 0-ZZnONI 
structure respectively contains more Zn and O atoms, and thus are unsaturated. However, 
spin-splitting effect can only be observed for the 0-ZZnONI structure as shown in Fig. 
2(c). And the calculated total spin {Stot) is 2. From the calculated energy levels shown in 
Fig. 2(c), we see that the energy difference between the highest occupied molecular orbital 
(HOMO) of spin-up and spin-down electrons {Ed (HOMO)) is only 0.13 eV, which is very 
small and indicates that the spin-splitting effect is very weak. The smaller energy gap in 
individual gaps of spin-up and -down electronic states is defined as E*, which is 0.18 eV 
and zero for the Zn- and 0-ZZnONI structure respectively. The zero energy gap indicates 
that the 0-ZZnONI structure is metallic at its ground state. And the small energy gap of 
0.18 eV indicates that the Zn-ZZnONI structure is chemically not very stable. However, the 
Zn- and 0-ZZnONI structures still might be fabricated using some special methods, since 
unpolarized hexagonal ZnO monolayers have been successfully produced already Gj. 

The deformation electron density and spin density are then further analyzed to study the 
origin of the magnetism in the 0-ZZnONI structure. From the deformation electron density 
shown in Fig. 3(a), we can see that all the O atoms get some electrons from Zn atoms in 
the nanoisland. Mulliken charge density analysis further shows that the O atoms at the 
corners and edges get about 0.1 e more electrons than the O atoms in the middle of the 
nanoisland. These extra electrons might supply the excessive electrons in one spin state. 
Detailed wavefunction analysis proves that the electronic states around the Fermi energy 
are mainly contributed by the oxygen-dominated edge states. Figure 3(b) shows the spin 
density distribution in the 0-ZZnONI structure. In this case, it is clear that the magnetic 
behavior is due to the oxygen edge states. 

The fact that the 0-ZZnONI structure without any magnetic impurities is magnetic might 
help us better understand the profound origin of magnetism in ZnO-based nanostructures. 
Besides, the magnetism in the 0-ZZnONI structure also hints that it might has direct appli- 
cations in nanoscale spintronics. So we further investigate the size effects on the magnetism 
in different 0-ZZnONI structures, which can be differentiated by the number of O atoms at 
each edge (z), and we will use Oi-ZZnONI to represent the 0-ZZnONI structure of different 
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sizes. Based on this definition, the previously discussed 0-ZZnONI structure is actually the 
06-ZZnONI structure. 

Table I shows the calculated magnetic and electronic properties for some Oi-ZZnONI 
structures. We find that the value of all (HOMO) ranges between and 0.2 eV, which 
is very small and indicates that the spin-splitting effect is weak for all the Oi-ZZnONI 
structures. From the calculated 5'tot listed in Table I, we can see obvious quantum size effects. 
Except for 03-ZZnONI, the Oi-ZZnONI structure with a larger size always has a larger 
total spin. Prom careful wavefunction analysis, we find that in the 03-ZZnONI structure, 
electrons of the middle oxygen atoms also contribute to its magnetic states. Considering 
that the magnetic states in other Oz-ZZnONI structures are all contributed by their O edge 
states, this explains why the 03-ZZnONl structure has a larger 5'tot than the 04-ZZnONI 
and 05-ZZnONI structures. However, the abnormal large total spin of the 03-ZZnONI 
structure is because of its too small size. For most Oi-ZZnONI structures, larger ones will 
have larger total spin. 

In summary, we have systematically investigated the structures and magnetic properties 
for 2D ZnO nanoislands. It is found that the structural reconstructions are mainly around 
the three corners and three edges for triangular ZnO nanoislands. Among the five ZnO 
nanoislands with different edge structures, only the 0-ZZnONI structure is magnetic at its 
ground state. And its spin density mainly distributes at the edge oxygen atoms. For the O- 
ZZnONI structure with different sizes, we reveal that the magnetic and electronic properties 
are quite dependent on their sizes. The Oi-ZZnONI structure with a larger size always has 
a larger magnetic moment, implying their potential usages in spintronics. 

This work was supported by the NSFC under grants No. 10604010 and No. 60776063. 
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TABLE I: The energy difference of HOMO (Ed (HOMO)), the smaller energy gap in individual 
gaps of spin-up and spin-down states (E*), and the total spin (^tot) for the Oi-ZZnONI structure 

with different sizes. 

i Ed (HOMO) (eV) E; (eV) S^ot 



3 0.01 0.19 2 

4 0.08 0.06 1 

5 0.17 0.06 1 

6 0.13 0.00 2 

7 0.01 0.02 3 
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List of captions 



Fig.l (Color online). Atomic structures from top view for the AZnONI (a), Zn-ZZnONI 
(b), 0-ZZnONI (c), BAZnONI (d) and BZZnONI (f). (e) and (g) Atomic structures from 
side view for the selected atoms in (d) and (f). In all the structures, red and grey balls 
represent O and Zn atoms respectively. 

Fig.2 (Color onhne). Energy levels of AZnONI (a), Zn-ZZnONI (b), 0-ZZnONI (c), 
BAZnONI (d) and BZZnONI (e). The left and right columns in each figure correspond to 
spin-up and spin-down electronic states, respectively. 

Fig.3 (Color online). Contour plot of the deformation electron density (a) and isosurface 
plot of the spin density (b) for the 0-ZZnONI. Red and grey balls respectively represent O 
and Zn atoms. The isosurface with the value of 0.01 e/ A? is shown in blue. 
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